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Abstract 

We show that the Density of States (DoS) for lattice Self Avoiding Walks can be estimated by using 
an inverse algorithm, called flatlGW, whose step-growth rules are dynamically adjusted by requiring the 
energy histogram to be locally flat. Here, the (attractive) energy associated with a configuration is taken to 
be proportional to the number of non-bonded nearest neighbor pairs (contacts). The energy histogram is 
able to explicitly direct the growth of a walk because the step-growth rule of the Interacting Growth Walk 
1^ samples the available nearest neighbor sites according to the number of contacts they would make. We 
have obtained the complex Fisher zeros corresponding to the DoS, estimated for square lattice walks of 
various lengths, and located the 9 temperature by extrapolating the finite size values of the real zeros to 
their asymptotic value, ~ 1.49 (reasonably close to the known value, ~ 1.50 |i50). 

Keywords: self avoiding walk; interacting growth walk; density of states; zeros of partition function; energy his- 
togram; theta point; monte carlo simulation 
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I. INTRODUCTION 



One of the challenging computational problems in the study of Self- Avoiding Walks (SAW) yj] 
- a statistical mechanical model for describing the excluded volume interaction of a linear polymer 
chain - is to estimate the configurational Density of States (DoS) as accurately as possible so 
that their equilibrium thermodynamic properties can be worked out. Equivalently, the problem 
is to find an efficient method for generating all possible configurations of a SAW. Since the total 
number of SAW configurations increases exponentially with the length, N, of the walk, there is a 
need to devise efficient sampling techniques for generating as wide a spectrum of configurations 
as possible. 

The standard unbiased methods of generating SAW configurations 121], on a regular lattice, not 
only suffer from sample-attrition problem but also fail to generate low-entropy configurations in 
sufficient numbers. In other words, an unbiased method generates only typical SAW configura- 
tions with equal a priori probability; low-entropy (or compact) configurations are rarely generated 
because of their very low probability of occurrence in the microcanonical ensemble of all possible 
configurations. So, Kinetic Growth Walk (KGW) yj] was proposed as a method for introducing 
local bias in the generation process so that the typical SAW configurations generated are more 
compact than the ones generated by the unbiased method; the sample-attrition problem also is less 
severe. The probability of generating a given KGW configuration depends on the way it is grown 
step-by- step. Since the cumulative bias introduced for generating a successful KGW configuration 
is known, it can be removed by suitably weighting the configuration. 

Recently, Prellberg and Krawczyk have shown that the KGW algorithm is actually an ap- 
proximate counting algorithm in which the weights associated with the walk configurations are 
the estimates of the number of configurations. During the generation process, decision to stop or 
continue the walk (a certain number of times) is based on how the weight, a^, associated with 
it compares with the average of weights, d^, associated with all the successful, as well as failed, 
configurations thus far generated. After a sufficiently large number of trials, gives a fairly 
accurate estimate of the microcanonical partition function for SAW. This algorithm, known as 
'flatPERM', can be used for estimating the DoS if the configurations are binned with respect to 
the number of contacts (i.e., non-bonded nearest neighbor pairs) they make. A reasonably flat 
contacts-histogram, over the full range of contacts possible for a given value of N, provides an 
independent numerical evidence that rare configurations are not missed. 
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An interesting inverse problem is to find out whether the requirement of a flat energy-histogram 
(or equivalently, contacts-histogram) could lead to an accurate estimate of the DoS. We show that 
such an inverse method can be devised and used for obtaining the DoS, taking the Interacting 
Growth Walk (IGW) algorithm as an example. IGW [|61] is a finite temperature generalization 
of the KGW that grows by Boltzmann- sampling the locally available sites according to the number 
of contacts (or, 'energy') they would make. We show that the inverse method dictates adjusting 
the parameter, which is equivalent to adjusting the weight of the growing configuration. We 
compute the Fisher zeros llsl-llin corresponding to the DoS, estimated for SAW on a square lattice, 
and locate the ^-point. Agreement with the expected value of the transition temperature justifies 
the idea of an inverse algorithm which is similar in spirit to the Reverse Monte Carlo (RMC) 
algorithm of McGreevy and Pusztai |17|]. 

The paper is organized as follows. In the next section, we briefly review the IGW model for 
generating lattice SAWs whose compactness can be tuned by a parameter. We point out that 
averaging the bias-correcting weights over all attempts to generate walks of given length gives an 
estimate of the microcanonical partition function. In the first part of section III, using the examples 
of KGW and IGW on a honeycomb lattice, we show that the growth rule of flatPERM consists of 
two independent parts - (i) a local KGW rule and (ii) a global weight-based rule for continuing 
a walk; in particular, we show that the weight-based rule cannot in general be reinterpreted in 
terms of a 'bath' temperature. In the second part of this section, we define the inverse algorithm, 
flatlGW, whose step-growth rules are dynamically adjusted by the requirement of local flatness 
for the energy histogram. In section IV, we present the Monte Carlo estimates of the DoS obtained 
for SAWs on a square as well as triangular lattice. We also present an extrapolation analysis of the 
Fisher zeros obtained for square lattice SAWs. In the last section, we summarize our results. 



II. INTERACTING GROWTH WALK (IGW) 

Starting from an arbitrary site of a regular lattice of coordination number, z, the first step is 
taken in one of the available z directions; choice for the subsequent steps is made from among 
the z — 1 nearest neighbor sites in the forward direction. In the unbiased sampling method, an 
attempt to generate a walk is discarded if the site chosen has been visited already. Even though 
this method ensures that the successful configurations are all generated with equal probability, 
pSAw _ _ -^^N^ (-j^g rejection rate increases exponentially with the length of the walk because 
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of which it fails to generate sufficient number of successful walks needed for good statistical 
analysis. 

An immediate solution to this problem is to choose at random one of the available, or unvisited, 
nearest neighbors for the next step to be taken. The number of available sites for the N^''' step, 
aAr(Civ_i) (e [0, z — 1]), depends on the configuration C^^i already generated. It is clear that 
this method, known as the KGW, introduces bias at every step of the walk. However, sample- 
attrition is less severe in this case and is only due to the geometrical 'trapping' {aN{CN~i) = 0) 
of the walk. The cumulative bias introduced during the generation process can be removed by 
associating a weight, A^^^{Cm) = z 11x^=2 '^k{C-k-i), with a configuration. These weights, 
referred to as the Rosenbluth-Rosenbluth (RR) weights [|3|], may be used for estimating the mean 
squared end-to-end distance. 

As becomes larger, the distribution of these RR- weights becomes flatter and so, the weighted 
averages fail to provide bias-corrected estimates in the asymptotic limit of very long walks. 

Averaging the RR-weights over all the configurations, trapped or successful, we obtain an 
estimate of the total number of A^-step SAWs 

ySAW 1 \^ aKGW/p \ /o\ 

where S is the total number of attempts made by the KGW algorithm to generate A^-step SAWs; 
since the weights associated with those configurations that have been trapped at various stages 
are zero, the summation is only over all the successful configurations. Zf/^^ will give the total 
number (equivalently, the microcanonical partition function) of A^-step SAWs, Zff""^ , in the limit 
iS — > oo. 

Since any given KGW configuration, Cn, will have a certain number of non-bonded nearest 
neighbor pairs {contacts), we can equivalently define. 



2^AW _ ap^^rn] 0'N,m = ^ ^n'^^ N , (3) 

m=l Cjv 

where A^'^^{Cn, rn) is the RR- weight associated with a configuration that has m contacts. In the 
limit 5 — 7- oo, therefore, aN^m is the total number of A^-step SAWs with m contacts, which is the 
DoS. It must be noted that the KGW algorithm is blind to the contacts being made during growth 
- i.e., identifying and counting the contacts are not a part of the basic growth rule. 



On the other hand, IGW algorithm |^ samples the available nearest neighbors on the basis of 
the number of contacts they would make. The probability that an available nearest neighbor site, 
say r^il), for the A^*'' step will be chosen is given by 

-el3Gm{rN{l)) ajv(Civ_l) 

P^^^(C^_i,/3g) = ^7^ a^(C^_i,/3G)= Yl e^'^'^ '"(^-W) (4) 

where m{rj^{l)) is the number of contacts the site rj^{l) would make, and e is the 'energy' cost of 
a contact. The parameter (3g is an inverse 'temperature' that helps sample the local environment. 
For (3g = 0, this is just the KGW growth rule; for (3g oo, only the site making a maximum 
number of contacts will be chosen. Interestingly, there is no sample-attrition in the limit /3g — ?• oo 
because a contocf-making site cannot be avoided during growth, and the walks thus generated are 
maximally (not globally) compact. Even though (3q^ cannot be identified with a canonical 'bath' 
temperature |113i|, it has been shown |^ that a coil-to-globule transition can be identified for some 
lattice-dependent value of 

As in the case of KGW, a weight Aj^^{CN, (5g) can be associated with a given IGW configu- 
ration, 

N 

I3g) = zx\[ a^(C^_i, /3g) (5) 

K=2 

SO that an estimate of the microcanonical partition function for SAW is given by the average. 



Cn 



where S{(3g) is the total number of attempts made by the IGW algorithm to generate A^-step 
SAWs, and miCN) is the number of contacts in Cn- In the limit S{I3g) — ?■ oo, the microcanonical 
partition function obtained this way is independent of the parameter I3g- The DoS is then given by 

aN,m = Jim Yl ^N^^i^N, /3g, m)e^'^« "^('^-) (7) 

<S(/3g)^oo <b{PG) ^ 

where A^j^"^ [Cn , I^g^^) is the weight associated with an IGW configuration Cn with m con- 
tacts. For (3g = 0, it is just the RR-weight. On the other hand, averaging only the weights, 

^\Pg) ~r 

Cjv 

leads to the canonical partition function at the 'bath' temperature in the limit S{(3g) — oo. 
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III. flatlGW - AN INVERSE ALGORITHM 



As mentioned earlier, a typical KGW configuration is more compact than a typical SAW con- 
figuration because it is generated by sampling only sites available for growth with the same prob- 
ability; this introduces a local bias not because it distinguishes between the available sites but 
because the number of available sites is random and configuration-dependent. flatPERM uses the 
RR-weights to decide whether the walk should be continued and with what adjustment to its RR- 
weight. So, the growth rule consists of two independent parts - (i) local (equivalently, the KGW) 
sampling of the available nearest neighbors with equal probability, and (ii) global check on how 
the RR-weight of the currently growing configuration compares with a cumulative average. The 
local (KGW) growth rule is independent of the global weights-based rule and, is never modified. 
The objective of this section is to use the solvable models of KGW and IGW on a honeycomb 
lattice in order to demonstrate that this independence is mutual. 



A. flatPERM versus step-growth rule 



(1) KGW: 



Poole et al, [|12l1 have shown that KGW on a honeycomb lattice is equivalent to the Interacting 
Self-Avoiding Walk (ISAW) at the inverse 'bath' temperature (3K{hc) = In 2. This is so because, 
on a honeycomb lattice, a walk makes a contact if only one nearest neighbor site is available for 
further growth. Therefore, the probability of generating a walk. Cat, with m{C) contacts is given 
by 

PI^^^{Cn) = - X x2"' = P^^^ X e-^'"^^^) '^^C'^) (9) 

3 2 

where each contact is assigned energy e = —1. It is clear that KGW is equivalent to ISAW at a 
'bath' temperature = (ln2)^^ Since the 'bath' temperature is higher than the temperature 
on a honeycomb lattice (~ 1), KGW is in the same universality class as SAW and therefore, 
unweighted averaging of the end-to-end distances of KGW also yields the SAW value for the size- 
exponent = 3/4. However, configuration-dependence of P^'^^{Cn) necessitates weighting the 
configurations appropriately so that they represent a bona fide statistical ensemble. 

The (microcanonical) weight associated with a KGW configuration. Cat, is then given by 

aT\Cn) = (P^^^(C^))-^ = '-^sAw- (10) 
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This is just the RR-weight for a KGW configuration because (3K{hc) = In 2 corresponds to the 
athermal case on a honeycomb lattice. The flatPERM parameter, r^, is then defined as the ratio 

,K{hc), 



i^^) AK{hc) _ ST^ 4^{/ic)/ 



- -,K{hc) ' =^Z^^JV \}-n) (11) 



where S is the total number of attempts made (successful or failed). Since P^^^ is a constant, we 
have the ratio, 

'^Af C'jv 

If ttv < 1, the walk will be continued with probability after resetting its weight to be equal 
to This situation is realized when the current configuration is more compact (with e = 

— 1) than an average or typical configuration. On the other hand, if r^v > 1, it is copied c = 
min([r7v], -2a) times, and each copy continued after reducing its weight by c; here, [r^r] is the 
largest integer less than and Za{= 1,2) is the number of available sites. In this case, the walk 
is less compact (with e = —1) than a typical configuration and therefore can be reused; 

Resetting the weight, a^'^^'^^ a^'^'^^'' = a^'^^'^^ or a^^^^'^ jc, is formally equivalent to resetting 
the inverse temperature PiK{hc) to a value given by 

?K[hc) P'Kihc) = l^r^lna^'-'^''^ (13) 

em[CN) 

This only provides an alternative interpretation of the global weight-based rule for a chain-growth 
process because of the known fact that the 'bath' temperature decides the degree of compactness 
of a typical configuration. However, this cannot be used as a parameter in the local growth rule 
because there is nothing to distinguish one available nearest neighbor from another for taking a 
given step. 

On a lattice of coordination number z > 3, even a configuration-independent bath temperature 

for KGW cannot be defined 1131] . For example , on a square lattice {z = 4), a walk makes a single 

(or double) contact if the site chosen has one (or two) occupied non-bonded nearest neighbor. 

Note that the walk is geometrically trapped if it chooses a site having three occupied non-bonded 

nearest neighbors. It is clear that the number of single (mi) and double (1712) contacts made is 

configuration-dependent. Averaging over all the configurations having the same total number (m) 

of contacts, but with different sets of single and double contacts, {(mi, 1122) | m = mi + 1712}, we 

can define an m-dependent inverse temperature, 

miln(3/2) + 2m,ln3 
m 
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where mi^2 denotes configurational averages. Since these /3-values are 



sharply distributed about 



the value f3sq ~ 0.45 (width ~ 0.02) in the asymptotic limit (A^ — )■ oo) UM, KGW may be said to 
be equivalent to ISAW at an ejfective bath temperature, /3sq ~ 0.45; but the converse is not true. 
If, however, the unweighted averaging of the end-to-end distances for KGW also gives the SAW 
value for z/(= 3/4) [|3D, it is because is higher than the 9 temperature (~ 1.15). Hence, even 
reinterpreting the global weight-based rule in terms of a 'bath' temperature is not straightforward 
for walks on a lattice of coordination number z > 3. 

(2) IGW: 

Since an IGW samples the available nearest neighbors on the basis of the number of contacts 
they would make, the individual steps of the walk on a honeycomb lattice may be identified either 



as two-indexed steps or as forced steps lll3h . A two-indexed step, labeled [mi, 1712], is a step to an 
available site making mi contacts rather than to the other available site making m2 contacts; on a 
honeycomb lattice, mi, m2 = 0, 1, 2. Steps labeled [2, 0], [2, 1] and [2, 2] indicate that the walk is 
geometrically trapped. Forced steps are those that have no other choice but to take the only site 
available; they are preceded either by contact-making two-indexed steps or by forced steps. 
The probability of generating a walk, Cn, with m{CN) contacts is given by 

"[0,l]+'^[l,2]+'^[l,0] / 9 \ "I0>2] 



where f3c is the inverse growth temperature that helps choosing an available nearest neighbor that 
makes a certain number of contacts; nj^^ ^2] is the number of steps that have chosen to make mi 
rather than m2 contacts. It is quite possible that different configurations with the same number 
of contacts m may have different sets of numbers {nymi,m,2\}^ we may use their average values, 
{^^[mi.ma]} ^ud define an inverse temperature, 

R (R m\ - R I ^[1-0] R (^[0.1] +^[1,2] +^[1.0]; , 

Pi{hc){PG,^> = PKihc)-\ TTT^Pg TTT^; in 



We may now rewrite Eq.lfTSl) in the form. 



m(C7v) em(Cjv) \l + e-'^o 

2^ 



mr^^i (16) 



For a given value of and (3g, the distribution of the inverse temperature l3i{hc)iPG,T^) over a 
large collection of IGW configurations is sharply peaked around an average value, say l3i(hc){l3G), 
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that may be taken to be the bath temperature. Unlike in the case of KGW, an implicit, 
configuration-independent, bath temperature does not exist for IGW > 0) on a honeycomb 
lattice. 

In order to remove the bias introduced in the walk-generation process, we associate the weight 
^n^^\Pgi Cn) with configuration Cn, 

AT\(iGM = (18) 

and compute the flatPERM parameter, r^, defined as 

3em(Civ) /3/{hc)(/3G: m{Cjv)) 

Again, depending whether < or > 1, the walk will be continued with a probability equal 
to rj\f or unity respectively; the weight a^'^''\f3G) will be correspondingly reset to the value 



a 



(Pg) = (/3g) or "tv Wg)/c where c = min(2, [tn])- Formally, this is equivalent 



N [Pg) — o-N \Pg) or 
to resetting fii{hc){PG, ^{Cn)) to a value, 

P'j^^^PgMCn)) = -^^\naT''\PG) (20) 

Since the partitioning of into the number of steps, ^j], that have chosen to make mi rather 
than m2 contacts (A^ = rijo^i] + ra^i^o] + ?t^[i,2] on a honeycomb lattice) is not unique for a config- 
uration Cat having a total of m (Cat) contacts, reinterpreting the global weight-based rule in terms 
of a 'bath' temperature is not possible even on a honeycomb lattice. 

However, by adjusting the value of (3g continuously during the growth process, we may be able 
to generate the desired SAW ensemble because, the local IGW rule distinguishes the available 
nearest neighbor sites on the basis of the number of contacts they would make. What is required 
is a global criterion that dictates appropriate adjustments to the value of (3g at every step. The 
fact that a flat energy-histogram is an indicator of an efficient DoS algorithm raises the question 
whether a requirement of local flatness for an energy-histogram, which is building up during the 
simulation, can be used as a global criterion for tuning /3g- 



B. flatlGW 



Let HximK) denote the number of A'-step IGWs having contacts {niK = 0, 1, 2, 



Ax) 

''K 



(x) 

where m)^' is the maximum number of contacts possible for a K-step walk). A perfectly flat 



9 



energy histogram corresponds to having the same value ifi^(m/^) for all uik. If HKimx) 7^ 
HK{mK — 1), the number of K-step walks having tuk contacts is more, or less, than those having 
{rriK — 1) contacts. So, the difference, 5HK{mK) = HK{mK) — HxirriK — 1), rnay be taken as 
a measure of the deviation from local flatness. 

Suppose, a Ji-step walk segment has already made mx contacts. Since it is an irreversibly 
growing walk, the number of contacts, m^-, cannot decrease during further growth; it can either 
increase or remain the same. Before taking the {K + 1)*'^ step, we should check whether it is 
desirable, or not, to increase the number of contacts. We assume, for simplicity, that no available 
site increases mx by more than one. 

Since /3c is the parameter that controls the extent to which a contact-making site is preferred, 
or avoided, its value for the {K + ly^ step may be chosen so that | 5HK+i{mK + 1) | does 
not increase. For example, if 5HK+i{rnK) > 0, (3g{K + 1) must have a negative value so that 
the local IGW growth-rule avoids a contact-making site; on the other hand, if 5HK+i{mK) < 0, 
I3g{K + 1) must have a positive value so that a contact-making site is preferred. The actual 
value of (3g{K + 1) should depend on | 5HK+i{mK + 1) |. The simplest choice is to have 
I I3g{K + 1) I oc I 5HK+i{rnK + 1) |, with the proportionality factor, say ^, chosen conveniently 
by trial and error. We have taken ,^ to be a uniform random number in the unit interval [14] and 
tried the following rule: 

I 1= min{| SHk+Mk + 1)^ I, I /^r^ 1} (21) 

where | | is the maximum value conveniently chosen to be within the overflow/underflow 

limits. The actual value of used for the {K + 1)*'^ step is given by 



Pg 



h/ A < 



- I Pg h/ a > where A = 5HK+i{mK + 1) ^ (22) 
z/ A = 0. 

Hence, the step-growth probability defined in Eq.© for constant [5g has to be redefined as 

Pk+1 ^^^^ K^^) 

where tk+i is the site taken by the {K + 1)*^ step, m{fK+i) is the number of contacts made by 
occupying that site (assumed to be either or 1 for simplicity) and a^+i is the number of sites 
available for the K + 1*^ step. The procedure, presented above, can easily be rewritten for the 
case when there is also an available site that increases rriK by more than one. 
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Depending on whether a contact is made or not at the {K + 1)*'' step, HK+iijriK + 1) or 
HK+iimx) is incremented by one count. Simultaneously, we have a variable gK+i{rnK + 1) or 
(7/^+i(mx) representing the cumulative weights of all the {K + l)-step configurations generated. 
The weights, essentially defined by Eq.© and Eq.®, have to reflect the fact that the growth 
parameter I3g{K + 1) is tuned at every step. 

The probability of generating a configuration Ck = {n, ^2, ■ ■ ■ , rR'} in the flatlGW algorithm 
is given by 

Ph''''{CK) = f{vi'''' (24) 

L=l 

Then, the variable gKi'mx) is updated with the weight corresponding to Ck' 

K 

9K{mK) ^ 9K{mK) + [PK^^iCKT^ ; rriK = J2 ^(^"'^) (25) 

L=l 

We obtain an estimate, gKirnx), of the DoS by averaging gKirriK) over all the attempts made to 
generate iiT-step walks. 

rv. RESULTS AND DISCUSSION 

In Fig.©, we show the reduced histogram, hxim) = Hxim) /max{HK{m)}, for walks of 
lengths in the range K = 25 to 300 generated on a square lattice. The value of (3^'^^ has been 
set equal to unity so that the adjustable parameter /3^^ G [—1, 1]. The flat histograms seen are 
evidence that the local criterion used for adjusting /3^\ Eq.(l2T] and Eq. dllj) . does bring about 
global flatness to the histogram data. The observation. Fig.©, that estimates of the DoS for 
short walks agree well with those obtained by exact enumeration gives confidence that the inverse 
algorithm, flatlGW, is as good as flatPERM. DoS estimated for longer walks, K in the range 10 to 
120 are shown in Fig.©. Reduced histogram and the DoS estimates obtained, with = 0.8, 
for walks on a triangular lattice are shown in Fig.©. Initially, a few trial runs may be necessary 
to fix suitable values of for different lattices. 

Estimates of the DoS can be used for computing the specific heat of ISAW at any given bath 
temperature and, in general, for computing the phase diagram. ISAW is known to undergo a 
transition from random coil phase to ^-phase and finally to a compact phase at temperatures that 
depend on the lattice used. So, locating a transition temperature, say the ^-temperature, will be a 
test for the accuracy of DoS estimates. 
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An interesting technique for locating a phase transition is to compute Fisher zeros of the canon- 
ical partition function, and check the (real) temperature to which they would converge. Since the 
walks generated are all of finite lengths, it is necessary to extrapolate the estimated temperatures 
to the asymptotic limit. 

Zeros of the polynomial, Zjv, in the variable x{N) (= e^) with the flatlGW estimates of the DoS 
as the coefficients, obtained for square lattice walks of length = 24, are shown in Fig.® along 
with those obtained by exact enumeration of 24-step SAWs on a square lattice. The agreement is 
quite good. As the walk length increases, it is expected that there will be a crowding of zeros in 
the neighborhood of the positive real axis; Fig.® demonstrates this crowding of zeros for walks 
of lengths N = 25, 55, 95 and 125. Numerical estimates of Xc{N), the value at which a real zero 
is expected, obtained for various values of are plotted as function of 1 /N in Fig.©. A simple 
extrapolation of the data gives the asymptotic value Xc{N — )■ oo) = 1.96 ± 0.01. This implies that 
/3g ^ 0.67 which is in close agreement with the known value ~ 0.667 1^. 

These results indicate that flatlGW can also be used for studying the statistical and thermody- 
namical properties of ISAW on a regular or disordered lattice. Recognizing the fact that the set of 
numbers, {///^(m)}, evolve as a rough surface, may be thought of as bounding lines within 
which the fluctuations of the iJ-surface are confined. So, it is intuitively clear that the iZ-surface 
will asymptotically reach a stationary 'flat' state; the time required for reaching this state and the 
asymptotic 'flatness' will depend on (3q"-^. 

flatlGW is an inverse algorithm in the following sense. The growth-rule of the flatlGW consists 
of two parts - (i) local IGW rule for sampling the available nearest neighbor sites and, (ii) global 
rule requiring that the energy-histogram evolves in such a way as to maintain local flatness. The 
latter is used for adjusting the parameter, /3g, on which the former depends. So, these two parts 
of the growth-rule are inter-dependent. An indicator that this algorithm works is the (statistical) 
convergence of DoS towards their stationay values, as the weights assigned to the walks are ac- 
cumulated independently during the simulation runs. In the direct methods, on the contrary, an 
indicator that the DoS converges to their stationary values is flattening of the energy-histogram 
during simulation; moreover, the local and the global components of the growth rule are indepen- 
dent. 

flatlGW is an inverse algorithm similar in spirit to the Reverse Monte Carlo (RMC) algorithm 
of McGreevy and Pusztai [[zl], in which the model-parameters are adjusted so as to be consistent 
with the experimental data rather than by the standard Metropolis criterion. 
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However, the local flatness criterion, used in the step-growth rule, is not such a rigid con- 
straint that a perfectly flat energy-histogram may be expected asymptotically. For instance, when 
I 5HK+i{mK + 1) 1= for all the available sites, (5q^^ = and therefore one of the sites will be 
chosen at random with equal probability; this will spoil the local flatness. 

It must be mentioned here that the rule, Eq.(l2TI). used for adjusting in the algorithm is ar- 
bitrary and introduces one more parameter, (3q°'^, besides using a randomizing factor ^. This will 
affect the computational efficiency of the algorithm when compared to the other ones such as flat- 
PERM [|4|] and Multi-canonical chain growth algorithm lilSll . Nevertheless, the main objective of 
this work is to demonstrate that an inverse, histogram-driven algorithm (a la RMC) can be con- 
structed for computing the DoS of linear polymers. A detailed comparative study of this algorithm 
vis-a-vis others will be helpful. 



V. SUMMARY 



We have demonstrated that an inverse algorithm, flatlGW, that is similar in spirit to the Reverse 
Monte Carlo algorithm [Q] can be implemented for estimating the DoS of linear lattice poly- 
mers. This algorithm generates a SAW step-by-step and irreversibly by using the IGW growth 
rule, Eq.dH), for a biased sampling of the available nearest neighbor sites; whether a contact- 
making available site is preferred or not is decided by a parameter, /^g, whose value is decided 
by the requirement of local flatness for the evolving energy-histogram. The DoS converge to their 
stationary values, as the statistical weights associated with the walk configurations accumulate 
independently during the simulation runs. 

We have used this inverse algorithm to estimate the DoS for SAW on a square, as well as trian- 
gular, lattice. We have located the Q temperature for square lattice SAWs by analyzing the Fisher 
zeros corresponding to their DoS estimates. Agreement with the known value lends numerical 
support to this idea of an inverse algorithm. 

*Author for correspondence. 
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FIG. 1: Plot of reduced histogram versus energy En and length n = 25 to 300 (insteps of 25: right 

to left) for ISAW. Energy —En is just the number of contacts made by a walk of length n. These have been 
obtained from flat histogram IGW on a square lattice with /3^"^ = 1.0 and M = 10^ attempts. 
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FIG. 2: Semilog plot of estimated density of states g'^^{En) versus energy En and length n = 5 to 25 
(insteps of 5) for ISAW on a square lattices obtained from fiatlGW with = 1.0; number of attempts 
M = 10*. Energy —En is just the number of contacts made by a walk of length n. The Monte Carlo results 
{•k) matches very well with exact (+) density of states. ( Monte Carlo errors are very small; maximum of 
error does not exeed one percent and are within the symbol size on the scale). 
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FIG. 3: Semilog plot of estimated density of states g'ff*'{En) versus energy En and length n = 10 to 120 
(insteps of 10: right to left) for ISAW on a square lattice obtained from flatlGW with /3^"^ = 1.0; number 
of attempts M = 10^. Energy —En is just the number of contacts made by a walk of length n. 
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(b) 

FIG. 4: Flat histogram IGW on a triangular lattice, (a) Plot of reduced histogram (b) Semilog plot of 
gn^{En) vs En for various n; Pq"'^ = 0.8; number of attempts Af = 10®. Energy —En is just the number 
of contacts made by a walk of length n. 



16 



N =24 





9 


® IS 


& 



Exact 
X MC 


® 8 



.3 I ^ ^ 1 1 1 1 

-3-2-1 1 2 3 

Real 



FIG. 5: Distribution of the roots of partition function of ISAW for A'^ = 24 on a square lattice. The results 
obtained from exact enumeration (o) and flatlGW with P^"-^ = 1.0 (x). We have not marked two of the 
roots since they lie far away from origin. 
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FIG. 6: Distribution of the roots of partition function of ISAW on a square lattice for (a) A^ = 25 (b) A^ = 55 
(c) A^ = 95 (d) Af = 125. The results obtained from flatlGW with /^^''^ = 1.0. Note that some roots, lying 
too far away from the origin, have been omitted. 
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FIG. 7: Results obtained from flatlGW on a square lattice with (3^^^' = 1.0. An extrapolation of the real 
part of leading zero, Xc{N), for = 11 to 51 having the smallest magnitude imaginary component near 
positive real axis. The intercepts of the linear fit gives Xc{N oo) = 1.9586 ± 0.011. 
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